Heat capacity at the glass transition 
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A fundamental problem of glass transition is to explain the jump of heat capacity at the glass 
transition temperature T g without asserting the existence of a distinct solid glass phase. Similar 
problems are also common to other disordered systems, including spin glasses. We propose that if 
T g is defined as the temperature at which the liquid stops relaxing at the experimental time scale, 
the jump of heat capacity at T g follows as a necessary consequence due to the change of the liquid's 
elastic, vibrational and thermal properties. In this picture, we discuss time-dependent effects of 
glass transition, and identify three distinct regimes of relaxation. Our approach explains widely 
observed logarithmic increase of T g with the quench rate and correlation of the heat capacity jump 
with liquid fragility. 



INTRODUCTION 

When a transition takes place between two distinct 
phases, the change of heat capacity and other thermody- 
namic quantities is consistently understood in a theory 
of phase transitions [1]. Often a disordered system such 
as a liquid forms a similarly disordered solid glass with- 
out a transition into a different phase, yet heat capacity 
changes with a jump. The jump is considered a hallmark 
of glass transition, and defines glass transition tempera- 
ture T g . The heat capacity jump immediately presents 
a problem that is at the heart of glass transition 0, Q : 
how can the jump be understood if there is no distinct 
second phase? 

This problem remains unsolved and controversial. One 
set of theories rationalize the jump of heat capacity by 
invoking thermodynamics of phase transitions. An in- 
structive illustration is the ongoing discussion of a popu- 
lar theory of glass transition, the Adam-Gibbs theory [J]. 
The theory connects the change of heat capacity at T g 
to the configurational entropy which becomes zero below 
T g where a phase transition between a liquid and a glass 
takes place [H, Q . The Adam-Gibbs theory has been con- 
vincingly criticized for a number of important reasons Q . 
Chief of these, also present in other similar theories, is 
that it has not been possible to identify the second low- 
temperature phase (the glass phase) . To circumvent this 
problem, several theories have subsequently put forward 
the proposals about the non-conventional mechanisms of 
the phase transition and non-trivial descriptions of the 
second phase, while retaining the idea of a phase transi- 
tion of sort |2| • 

Another set of glass transition theories consider that 
glass transition phenomena at T g have purely dynamic 
origin, and simply correspond to the freezing of atomic 
jumps in a liquid at the experimental time scale Q. The 
absence of a phase transition and thermodynamic effects 
at T g is supported by the wide experimental observation 
that the liquid and the glass at T g have nearly identical 
structure [2]. However, the challenge for the dynamic 



theories is to explain both the origin of the jump of heat 
capacity at T g and its large magnitude, which for some 
systems can be of the order of &b per atom. 

In addition to glass transition in structural liquids, 
similar problems exist in other disordered media. For ex- 
ample, spin glasses have seen large developments of ideas 
based on phase transitions and the existence of the sec- 
ond distinct phase. Similar to the structural glass transi- 
tion, these theories have been used to explain the cusp in 
susceptibility at the glass transition temperature. Sim- 
ilar to the structural glass transition, several important 
problems remain in this field as well, including identify- 
ing the nature of a distinct spin glass phase, dependence 
of the cusp on field frequency or observation time, slow 
relaxation effects etc @. 

In this paper, we propose how to explain the jump of 
heat capacity in a purely dynamic picture, without as- 
serting the existence of a distinct thermodynamic solid 
glass phase and, therefore, show how to reconcile the 
above contradiction. We recall that glass transition tem- 
perature T g has two experimental definitions which give 
similar values of T g . In the calorimetry experiments, T g is 
the temperature at which the jump of constant-pressure 
heat capacity, C p , is seen. In the experiments that mea- 
sure t (e.g. dielectric relaxation experiments), T g is the 
temperature at which r exceeds the time of the exper- 
iment t of about 10 2 — 10 3 s. We propose that when 
r exceeds t, the jump of heat capacity at T g follows as 
a necessary consequence because freezing of local relax- 
ation events alters liquid elastic, vibrational and thermal 
properties including bulk modulus and thermal expan- 
sion. Hence, there is no need to invoke the existence of a 
second glass phase and a phase transition of sort. In this 
picture, we discuss time-dependent effects of glass tran- 
sition, and identify three distinct regimes of relaxation. 
Our approach explains widely observed logarithmic in- 
crease of T g with the quench rate and the correlation of 
heat capacity jump with liquid fragility. 
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CHANGE OF HEAT CAPACITY AT T g 

The commonly discussed quantity from the calorime- 
try experiments is the ratio of constant-pressure liquid 
heat capacity, CL to glass heat capacity, Ci. We do 
not consider the overshoot of heat capacity on heating 
and its undershoot on cooling, discussed elsewhere @, HI ■ 
Heat capacities are considered at temperatures separated 
by the interval in which these effects decay to the val- 
ues of Cp and C| attributed to the liquid and the glass 

[3,0]. For various liquids, §| = 1.1 - 1.8 [3]. Using 
the known relationship C p — C v = VTa 2 B, where C v is 
the constant- volume heat capacity, a is the coefficient of 
thermal expansion and B is bulk modulus, we write 



ViT l afB l 



V g T g o? g B g 



(1) 



where i and g subscripts refer to the liquid and the glass. 

We start with addressing the origin of the difference 
between Bi and B g and between a; and a g . Unlike in 
a solid glass, atoms in a liquid are not fixed, but rear- 
range in space. This gives liquid flow. Each flow event is 
a jump of an atom from its surrounding cage, accompa- 
nied by large-scale rearrangement of the cage atoms. We 
call this process a local relaxation event (LRE). A LRE 
lasts on the order of Debye vibration period r = 0.1 
ps. Frenkel introduced liquid relaxation time r as the 
time between LREs at one point in space in a liquid [l(| ■ 
Frenkel's main idea was that at short times t < r, liq- 
uid response is the same as that of a solid, i.e. is purely 
elastic. On the other hand, for t > r, viscous flow takes 
place, during which each LRE is accompanied by addi- 
tional, viscous, displacement. Hence, for t > r, liquid 
response to external perturbation (e.g . pressure) con- 
sists of elastic and viscous response [lfj . This discussion 
provided the microscopic basis for the earlier phenomeno- 
logical model by Maxwell [ll[ > who proposed to separate 
elastic and viscous response in his viscoelastic approach 
to liquid flow. 

Lets consider that pressure P is applied to a liquid. 
Pressure induces a certain finite number of LREs, which 
bring the liquid to the equilibrium state at new external 
conditions (P, T) after time r. Following the Maxwell- 
Frenkel approach, the change of liquid volume, v, is 
v = v e i + v r , where v e i and v r are associated with solid- 
like elastic deformation and viscous relaxation process 
due to LREs, respectively. Lets now define T g as the 
temperature at which r exceeds the observation time t. 
This implies that LREs are not operative at T g during the 
time of observation. Therefore, v at T g is given by purely 
elastic response as in elastic solid. Then, P = Bi v ^ l y^ r 

and P — BgY^ where VJ° and V g are initial volumes of 
the liquid and the glass and v g is the elastic deformation 



of the glass. Let AT be a small temperature interval that 
separates the liquid from the glass such that r in the liq- 
uid, t u is n = r(T g + AT) and < 1. Then, V° » V° '. 
Similarly, the difference between the elastic response of 
the liquid and the glass can be ignored for small AT, 
giving v e i w v g . Combining the two expressions for B\ 
and B g , we find: 



Bj 



B„ 



ei + 1 



(2) 



where ei — ^f- is the ratio of rclaxational and elastic 
response to pressure. 

u\ can be calculated in a similar way. Lets consider liq- 
uid relaxation in response to the increase of temperature 
by AT. We write at = yr Vc ^ r and a g — yw-^, where 
v e i and v r are temperature-induced volume increases in 
a liquid that are related to solid-like elastic and relax- 
ational response, respectively, and v g is the elastic re- 
sponse of the glass. Combining the two expressions and 
assuming VJ° = V g and v e i = v g as before, we find 



on = (£2 + l)a g 



(3) 



where e? = — is the ratio of rclaxational and elastic 

Z Vel 

response to temperature variation. 

Eqs. (|2|3j) describe the relationships between B and a 
in the liquid and the glass due to the presence of LREs in 
the liquid above T g and their absence in the glass at T g , 
insofar as T g is the temperature at which t < r. We note 
that relaxational response v r of both liquid and glass de- 
cays during time r. Because B g and a g correspond to 
t < t, B g and a g are unrelaxed, or non-equilibrium val- 
ues of bulk modulus and thermal expansion, respectively. 
This point is discussed in Chapter 3 in detail. 

Using Eqs. |03) in Eq. (JTJ), we write: 



a 1 



1 + ~/a g T g 



(4) 



where e = , 7 = Vg a gB g /CQ is the glass Griineisen 

parameter and C~* is the unrelaxed, or non-equilibrium 
(see Chapter 3) heat capacity of the glass. We note that 
in Eq. Q, as in Eqs. (03]), we set T t rj T g and V t w 

c' 

Vg. This underestimates the experimental setting its 
lower limit, because C p is measured in the finite range of 
temperature and volume such that T; > T g and Vi > V g . 



We now calculate 



Close to T g , C v is due to the 



vibrational motion only, whereas the contribution to C v 
due to the diffusional motion is negligible. This is an 
important assertion that perhaps was not appreciated 
before, and simplifies the problem greatly. The asser- 
tion follows from the explicit calculation of liquid C v as 
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a function of r [12| . and is consistent with the experimen- 
tal results showing that liquid C v is close to 3N around 



the melting point [13j, LLJ] . It also follows from a more 
general argument that does not rely on the explicit cal- 
culation of C l v , as discussed below. 

Above T g , each atom participates in the vibrational 
motion during time t and in the diffusional motion when 
it jumps between two equilibrium positions during time 
of about To- For t ^> r, atoms can be separated into two 
ensembles of N v n, vibrating and Ndif diffusing atoms. 
The energy of an atom in each ensemble consists of ki- 
netic and potential energy. Then, the partition sum 
is Z = Z V ib ■ Zdif, where Z v u, is due to vibrations 
and Z^f is due to diffusion. Liquid energy is E = 
T 2 Wr (^ n (Zvib • Zdif)) = T 2 -^p lnZ vib + T 2 -^ \nZ d if — 
E V ib + Edi f , where E v n, and Edi / is the energy of vibrat- 
ing and diffusing atoms, respectively. Here and below, 
temperature derivatives are taken at constant volume. 
At any given moment of time, Ndif = No exp(— U/T), 
where No is the total number of atoms in a system 
and U is the activation energy barrier for a LRE (U 
can be temperature-dependent). Combining it with r = 
To exp(t//T), where To is the Debye vibrational period of 
about 0.1 ps, we write 



N, 



dif 



No 



(5) 



which also directly follows by noting that the jump prob- 
ability is To/r. 

At T g , 231 10~ 16 , i.e. the number of diffusing atoms 
is negligible. Therefore, the ratio of the energy due to 
diffusion to the total energy, — is negligible. Similarly, 

< 1, giving 



dT 



hxZ, 



dT 



InZ, 



dif 

«; l 

vib 



(6) 



Liquid entropy is S 
(T-^p In Z vib + In Z vib ) 4 



= (Tln(Z V ib ■ Zdif)) 
(T-^p In Z dlf + In Zdif) 



^dT ^ n Zvib 



InZ, 



dif, 



In Z v ib 

have used Eq. @. Then, C' v = T^ = 



where we 

dS 



^dT V^dT m Zvib 

rp d (rp _d_ 1 7 

dT V dT ^vib 



^dT m Zvib 



T-^p\nZ vi b, where we have used 
Eq. © once more. Therefore, C v around T g (as well as 
at any temperature T such that <C 1) is essentially 
due to the vibrational contribution to Z. 

The vibrational states of a liquid are given by one lon- 
gitudinal mode and two transverse modes with frequency 
ui > 1/t [13]. If — -C 1, as is the case around T g , trans- 
verse modes in a liquid account for essentially all trans- 
verse modes that exist in a solid glass. Together with 
the fact that the phonon density of states increases as 
oc uj 2 , this means that the energy of the missing trans- 
verse waves with frequency u> < 1/t is negligible. Hence, 



the vibrational energy of a liquid in the regime - < 1 
can be calculated as the energy of all 3N phonons as in a 
solid glass 12]. Therefore, in discussing the vibrational 
C v of a liquid around T g , we can use the results derived 
for solids. 

The partition function of a harmonic solid is Z = 

(jF0 (^b=1), where uj is the geometrically aver- 
aged phonon frequency [l|, giving the free energy F — 
3NT\n^. In the purely harmonic case, to is constant, 



In J^- ) and 



giving the entropy S = — (ffr) = 3 AT (l 
C v = T (•§|r) = 3N . On the other hand, anharmonic- 
ity, particularly large in liquids, results in the decrease 
of uj with temperature. Importantly, as we show be- 
low, this decrease is different below and above T g be- 
cause a is different (see Eq. If uj is not constant, 



S = 3N (l + hi 



Cy = 3N 1 



T d 
u dT 



), and 



2T duj T 2 




2 T 2 d 2 uj 


uj dT uj 2 




uj dT 2 



(J) 



where the derivatives are taken at constant volume. 

The effect of anharmonicity can be discussed in 
the quasi-harmonic approximation by introducing the 
Griineisen parameter 7 = — ¥j (f^) T to the phonon pres- 
sure, P ph = - (|f ) T = ^P-. Then, the bulk modulus 
due to the (negative) phonon pressure is B p h 
and y 9 QT h J — where we neglected the depen- 

dence of 7 on V. Using 7 = V £ B and B = B + 
B p h, where B is the zero-temperature bulk modulus, 

9B ph 



3NTf 
V 



( dB ph \ 
V «T ) 



—a(B + Bph), where we set C v = 3N in 
this approximation. For small aT, which is often the 
case in the experimental temperature range, this implies 
B oc — T, consistent with the experiments [15j . We note 
that experimentally, B linearly decreases with T at both 
constant volume and constant pressure [15j . The de- 
crease of B with T at constant volume is due to the intrin- 
sic anharmonicity related to the softening of interatomic 
potential at large vibrational amplitudes; the decrease 
of B at constant pressure has an additional contribution 
from thermal expansion. Assuming uj 2 cx Bq + B p h and 

combining it with (J-§fr-^J — ~a(B + B p h) from above 
gives — — |r. Putting this in Eq. ([7]) gives 



C v = 3N(1 + aT) 



(8) 



We see that the derived expression for C v is linear with 
T and depends on a but not on uj, unlike in Eq. (|7J). This 
result follows from Eq. ([7]) as long as 4I? oc B or 4^ oc uj. 

From Eq. ©, §| = \Xl% at T g- Usin S [t in Ec l- © 
and retaining only linear terms in aT (we find that the 
linearization and direct combination of Eqs. (Q} and © 

C l 

give close values of below), we write 



4 



C l 

^L = l + iag T g (e-l)+T g ( ai -a g ) (9) 



Recalling that e = 



(£2+1) 



= 7? TT ' we see that Ec l- ® 



relates ^f- to the changes of a and B due to the presence 
of relaxational response in the liquid and its absence in 
the glass. 

Interestingly, Eq. ((9]) predicts that temperature de- 
pendence of C p should follow that of a. This is in agree- 
ment with recent simultaneous measurements of C p and 
a showing that both quantities closely follow each other 
across T g [l6(. 

Importantly, the jump of heat capacity at T g in our 
theory takes place within the same single thermodynamic 
liquid phase, but below and above T g the liquid has dif- 
ferent values of a and B due to the freezing of LREs at T g 
where the liquid falls out of equilibrium. In this sense, 
our theory is purely dynamic. In contrast to previous 
glass transition theories [2J , we do not discuss transitions 
between distinct thermodynamic phases, even though it 
may be tempting to invoke thermodynamic phase transi- 
tions, conventional or unconventional, in order to explain 
the heat capacity jump. 

c l 

In Table 1, we show for several common glass- 
es' 

formers with both small and large -^f in the range 1.1-1.8 



a g , ai, B g 

and Bi [l7H29| , we calculate ^ using Eq. ([9]) . Given the 



Using the experimental values of 7, T g 



approximations used, including T; — T g and Vi = V g that 

c l 

underestimate , Table 1 shows a reasonable agreement 
between the calculated and experimental values. The 
worse agreement for CKN is probably due to the fact 
that it is a solution Q, for which our approximations 
are expected to be less successful. We further remark 
that the agreement is subject to uncertainties of 7, a, B 
c l 1 — I 1 — I 

and 15, 21] which are, moreover, taken from different 
experiments. For these reasons, we view Table 1 as an 
illustration that the differences between the existing val- 
ues of a g and a\ and between B g and Bi are large enough 

C l 

to give the right magnitude of experimental 



general case of arbitrary relationship between r and t. 

c l 

This gives time-dependent properties of ^w- 

Lets now consider two liquids at two different temper- 
atures T\ and T2 with relaxation times T\ and t 2 such 
that T2 < T\ and t\ < T2, and calculate the ratio of 
their heat capacities, (f 1 ^- The response of both liq- 
uids to pressure now includes viscous relaxational com- 
ponent due to LREs. Hence, we write P = B\ 1 

and P = B 2 v " l y V o r ' 2 , where B ll B 2 are the bulk moduli 
and u r ,i, v r .2 are the relaxational responses of the two liq- 
uids, respectively. Combining the two expressions gives 



Bi 

B 2 



1 + 



"r,2 



1 



Vr.l 



(10) 



Similarly, considering temperature- induced response in 
the two liquids gives 



a 2 



1 



1 



(11) 



where ol\ and 012 are thermal expansion coefficients of the 
two liquids. 

When relaxational response in the low-temperature liq- 
uid is absent, u,., 2 = 0, Eqs. (llOHllj) become Eqs. ([2]|3]). 
Note that v r and v e i in Eqs. (|10M11[) are not the same 
because they are due to different effects of pressure and 
temperature. This difference will be accounted for below. 

We now recall that relaxation of liquids at low tem- 
perature follows slow stretched-exponential form: q = 



9o 1 



where q is a relaxing quantity, qq is its 



amplitude, t is observation time and /? is a stretching 
exponent [3, d, HO] • /3 decreases from 1 at high temper- 
ature to 0.5-0.8 at T g [3l|. Recently, we proposed that 
slow relaxation in liquids is a result of elastic interaction 
between LREs via high-frequency elastic waves that they 
induce [321 ] . The stretched-exponential relaxation follows 
as a result of this interaction [33] ] . 



Hence, 



= 9o 



1 



and 



V r ,2 



1 



where fi\ and /?2 are stretching ex- 



ponents in both liquids. Therefore, as follows from Eqs. 
(|9l lTT|) . time dependence of -^f^ is given by the following 
three equations: 



TIME-DEPENDENT EFFECTS 

The jump of heat capacity at T g in Eq. ([9]) is due to 
different a and B below and above T g due to the freezing 
of LREs at T g . This reflects the empirical definition of 
T g as the temperature at which r exceeds the observa- 
tion time t, as in a glass transition experiment. We now 
remove the empirical constraint r > t, and consider the 



= I+7CX2T2 I — — 1 



a 



p.2 



Bi 

B 2 



ol\ B 2 



1 + ei 1-e 



a i 

- tt2T 2 ( 1 

a 2 



1 + ei 1-e 



-(*) 



f>i 



(12) 



(13) 





7 


r 9 


a t ■ 10 4 


a 9 ■ 10 4 




S 9 




p 

T7W 

°P 


c 

p 

U P 






[K] 


[K- 1 ] 


[K- 1 ] 


[GPa] 


[GPa] 


calc. 


exp. 


Glycerol 


2.2 [17] 


190 [18] 


5 [19] 


1 [19] 


5.5 [20] 


9.9 [20] 


1.6 


1.8 [3] 


PVAC 


0.6 [21] 


304 [21] 


7.1 [21] 


2.8 [21] 


2 [21] 


3.5 [21] 


1.3 


1.4 [21] 


OTP 


1.2 [22] 


241 [18] 


7.2 [23] 


3 [23] 


2.1 [23] 


3.7 [23] 


1.3 


1.5 [3, 24] 


OTP-OPP 


1.3 [22] 


235 [25] 


8.5 [25] 


2.5[25J 


2.9 [25] 


5 [25] 


1.6 


1.5 [25] 


PS 


0.5 [25, 26] 


355 [25] 


6 [25] 


2.5 [25] 


1.5 [25] 


2 [25J 


1.2 


1.3 [25] 


CKN 


0.9 [21] 


340 [21] 


3.6 [21] 


1.2 [21] 


7.6 [21] 


15.9 [21] 


1.2 


1.6 [21] 


B2O3 


0.3 [21, 27] 


550 [21] 


4 [2Lj 


0.5 [21] 


2.6 [21] 


10 [21] 


1.3 


1.4 [21] 


NaAlSi 3 8 


0.35 [28] 


1100 [29] 0.54 [28] 


0.23 [28] 


20 [28] 


40 [28] 


1.05 


1.11 [28] 


Ge0 2 


0.27 [28] 


580 [28] 


0.76 [28] 


0.27 [28J 


8.08 [28] 


23.87 [28J 1.04 1.08 [3, 28] 



TABLE I: Experimental and calculated values of -^f. When more than one value of a, B, C v or 7 was available, we have 

taken the average. PVAC, PS, OTP, OPP and CKN stand for polyvinylacetate, polystyrene, orthoterphenyl, orthoterphenyl 
phenol and Ca(NOa)2-KN03 mixture, respectively. The ratio of the last term in Eq. ((9]), T g (ai — a g ), to the second term, 
■ya g T g {e - 1), varies from 0.1 in glycerol and 0.3 in OTP-OPP to 2 in NaAlSi 3 8 and 4 in GeQ 2 . 



I + C2 f 1-e ($)) 

_I = V L (14) 

- l + -."(*)*) ( } 

where we introduced relaxational amplitudes e± and e 2 
as in Eqs. 

There are three times in Eqs. (|12H14I) : t, t\ and r 2 . 
This sets three distinct regimes: (1) fast regime t <C t\ <C 
t 2 ; (2) intermediate glass transition regime Tj « t < t 2 
and (3) slow regime t\ <C t 2 <C <• Each regime sets a 
different mechanism governing the relationship between 
the heat capacities of the two liquids. 

Regimes (1) and (3) both give — = 1 and -g 1 = 1 in 
Eqs. (|13I14|) . and therefore C pA =C Pt2 from Eq. (TT21 . 
albeit for different physical reasons as discussed below. 

Regime (1) corresponds to non-equilibrium states of 
both liquids and to non-equilibrium values of their C p . 
In this regime, the terms in brackets in Eq. (fT3"|) and Eq. 
(fl4)) are close to because too little time elapses for any 
relaxation to take place. Physically, zero relaxational 
response in both liquids directly follows from Frenkel's 
idea that at short times, the response of a liquid is the 
same as that of a solid, i.e. purely elastic [10|. In this 
case, C P: \ — Cp.2 follows from the absence of relaxation 
in both liquids. 

Regime (3) corresponds to equilibrium states of both 
liquids and to equilibrium values of their C p . In this 
regime, the terms in brackets in Eq. (|13[) and Eq. (|T4|) 
are close to 1 due to relaxational response acquiring its 
maximal value as both liquids reach their equilibrium 
state, giving C Pi i = C Pi 2- Importantly, our theory pre- 
dicts no difference between C Pt \ and C Pt 2 in equilibrium 



(except for the residual difference due to T\ > T 2 and 
V\ > V2 in Eq. (JTJ) discussed earlier and the difference 
in C v in Eq. ([8]) related to a finite temperature interval 
where C Pi i — C Pj2 is measured). 

Regime (2) corresponds to the laboratory glass transi- 
tion when the first liquid is in equilibrium (ti <C t) but 
the second liquid is not (i <C T2). If the second liquid is 
defined as the glass, the change of heat capacity between 

-(j-Y 1 

the liquid and the glass follows: 1 — e v T i / is close to 

1 but 1 — e v^y i s close to in Eqs. (|13H14[) . giving 
B x = -^fj and a x = a 2 (e 2 + 1) as in Eqs. (gJE]) and 
therefore non-zero C Pi i — C P;2 in Eq. (fT2"j) . as in Eq. ©. 
This is the important result of our theory. 

We note that, by definition, the time range of regimes 
(l)-(3) is determined by the relationship between r and 
the observation time t. In the calorimetry experiments 
in particular, C p is measured on the time scale t = 
10 2 — 10 3 s. Depending on r, regimes (l)-(3) can be 
too fast or too slow for the experiment. For example, 
if ti,t 2 < 100 s, regimes (1) and (2) are too fast the 
experiment. On the other hand, regime (3) can be too 
slow for the experiment, and can last for astronomical 
times and longer. Lets consider Si0 2 glass at room tem- 
perature T r =300 K. The activation energy barrier U can 
be assumed temperature-independent, because Si0 2 is 
a strong liquid. Then, from r(T g ) — to exp(U/T g ) and 

Is. 

r(T r ) = T exp([//T r ), r(T r ) = r (^) . Taking 

t =0.1 ps, T g w 1500 K and r(T 9 ) = 10 3 s, r(T r ) = 10 67 
s, approximately the fourth power of the age of the Uni- 
verse. For t > r(T r ), Si0 2 at room temperature is an 
equilibrium liquid, and shows no jump of heat capacity 
on cooling from high temperature related to the freezing 
of LREs. 
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Further, our theory explains two widely observed and 
important effects of glass transition. First, a well-known 
effect is that T g , defined as the temperature at which the 
jump of heat capacity takes place in the calorimetry ex- 
periment, logarithmically increases with the quench rate 



q (see, e.g., Refs. [34j, |35[). According to our discussion 



above, the jump of heat capacity at T g takes place when 
the observation time t crosses liquid relaxation time r. 
This implies that because q = ^r- , r at which the jump 
of heat capacity takes place is r(T 9 ) = where AT is 
the temperature interval of glass transformation range. 
Combining this with r(T g ) = t exp(U/T g ) (here U is 
approximately constant because r is nearly Arrhcnius 
around T g 32]), we find: 



T a = 



U 



In AZ 

TO 



In q 



(15) 



According to Eq. (fl5|) . T g increases with the loga- 
rithm of the quench rate q. In particular, this increase is 
predicted to be faster than linear with In q. This is con- 
sistent with the experimental results [34j |. We note that 
this theory predicts no divergence of T g because the max- 
imal physically possible quench rate is set by the minimal 
internal time r , Debye vibration period, so that Ap in 
Eq. (fl~5j) is always larger than q. 

Second, Eqs. (|12H14[) predict the correlation of 
with liquid fragility @, HI]. Lets apply Eqs. (I12j|14j) 
to the glass transition regime where t <C r 2 . Then, 

1 — e V T 2/ = in Eqs. (|13M14|) . corresponding to the 
absence of relaxational response in the glass. Next, lets 
consider that AT g is the temperature interval in which 

c l 

£§■ is measured so that t\ — r(T g + AT g ) and t 2 « T~(T g ). 
Then, because /3 is nearly constant near T g [31], / = 



1-e" 
as / ? 



\ T1 ' in Eqs. (|12til4|) can be expanded around T g 
where the derivatives are taken around 

dlos 



T g . From the definition of fragility m = 



dr I 

dT\T=T g 
AT, 



r(T g ) 



to. Next, 



df I 



r{T g ) 



\T=T„ 



giving 



/ = -j^-m, and Eq. (|T2l becomes 



C l 

-4 = 1 + ja g T g 



^(^TO+1 

ei^TT-^TO + 1 



- 1 + a Q ATge 2 m 



If toA2: > 1, Eq. {T6j simplifies to 



(16) 



C p 
r<9 



i + M„Ar, ( . 2 (7^ + l] ,„ ( l.7i 



Eqs. (|16I17I) predict that increases with liquid 
fragility, provided other parameters do not significantly 



change. This is consistent with experimental data from 
a large set of liquids d, (3|| . For a wider range of chem- 
ically and structurally different liquids, the correlation 
holds within distinct families 37, 381. 



SUMMARY 

In summary, we observe that when a transition takes 
place between two distinct phases, anomalies in the ther- 
modynamic functions are explained by the phase tran- 
sition theory, a well-understood topic in physics If 
a distinct second phase can not be identified as in glass 
transition, the apparent anomalies can be explained in a 
picture that does not invoke phase transitions and ther- 
modynamics, but where the system stops relaxing on the 
experimental time scale. Our approach explains time- 
dependent effects of glass transition, including widely ob- 
served logarithmic increase of T g with the quench rate 
and the correlation of heat capacity jump with liquid 
fragility. 

We are grateful to M. Moore for discussions and to 
EPSRC and RFBR for support. 
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